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Abstract 

in 
in 

Previous results on the scattering of surface waves by vertical vorticity on 

(N ' 



shallow water are generalized to the case of dispersive water waves. Disper- 
sion effects are treated perturbatively around the shallow water limit, to first 
order in the ratio of depth to wavelength. The dislocation of the incident 



S3 

wavefront, analogous to the Aharonov-Bohm effect, is still observed. At short 

3 ' 

wavelengths the scattering is qualitatively similar to the nondispersive case. 
At moderate wavelengths, however, there are two markedly different scatter- 
ing regimes according to wether the capillary length is smaller or larger than 
\/3 times depth. The dislocation is characterized by a parameter that de- 
pends both on phase and group velocity. The validity range of the calculation 
is the same as in the shallow water case: wavelengths small compared to vor- 
tex radius, and low Mach number. The implications of these limitations are 

carefully considered. 
03.40.Kf, 47.35.+i, 47.10.+g 
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I. INTRODUCTION 



In a preceding paper |T[], hereafter refered to as I, we studied the scattering of surface 
waves by a stationary vertical vortex in the long wavelength approximation: surface tension 
was neglected and the fluid depth was supposed to be small compared to wavelength. This 
is also called the shallow water approximation. There were two motivations for the study 
of shallow water waves scattering. First, they are non dispersive waves, like acoustic waves 
in fluids, and it was plausible that a generalization of calculations for sound scattering by 
vorticity [|3| was feasible. Secondly, it was a first attempt towards a quantitative confirmation 
of the heuristic approach of Berry et al. 0. The aim of this paper is to go beyond this 
approximation. 

In actual experimental situations 0] the shallow water limit is hard to obtain and, if 
a quantitative comparison with experiment is desired, it becomes necessary to take into 
account the finite depth and the surface tension. The main difference between surface waves 
in shallow water and in deeper water lies in the fact that in the latter case dispersion 
effects are important: there are two length scales, one associated with depth and the other 
with surface tension, that are responsible for wave velocity depending on wavelength. In 
this paper we seek to describe the scattering of surface waves by vorticity in terms of a 
single differential equation in which surface elevation is the only dependent variable. This is 
possible in a perturbative treatement away from the shallow water case, and we here present 
results that correspond to first order corrections. 

As in I, we consider the scattering of surface waves by a stationary vortex, in the limit of 
small Mach number (the velocities of fluid particles are small by comparison with the phase 
velocity of the waves), M < 1, and large wavenumber k, i.e. (3 = ka ^> 1 where a is a 
typical length associated with the vortex flow. The product Mf3 is assumed to be of order 1. 
In Sec. [II], we derive, from the full hydrodynamic set of equations, an approximation valid 
to order O(M) (or 0{f3~ 1 )). First, equations are linearized for small surface perturbations 
around a steady vertical vortex and then higher order terms in M and (3~ l are discarded. 
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We shall pay particular attention to the orders of magnitude of the different terms, and 
will justify the neglect of dissipative effects. The recovery of the shallow water results is 
subtle since it involves taking the singular limit of vanishing surface tension. There appears 
a partial differential equation (Eq. ( |2.38| ) below) that contains a squared Laplacian, and it 
is reduced to our previous result (1-14) of ref. [I] in the shallow water limit, i.e. when the 
layer's depth is small and surface tension is negligible. 

The solution of equation Q2.38Q is given in Sec. |T|. The results, given by ( j3~7[) , (|3.27|) and 
( |3.28| ) and the calculations of the Appendix, seem much more complicated than the similar 
shallow water results (1-20), (1-24) and (1-25) of [I]. However, this complexity is essentially 
algebraic, and actually the physical results are rather similar, except when dispersive effects 
are closely balanced by advection to yield a spiral pattern for the scattered waves. The 
wavefront dislocation is characterized by a parameter a which is a generalization of the one 
in [I], and tends towards it smooothly in the shallow water limit. In the dispersive case, 
a depends on both the phase and group velocity of the waves. We give a perturbative 
justification of the heuristic argument of Berry et al. [^]. The behaviour of the scattered 
wave however depends strongly on the ratio of depth to capillary length. We also exhibit 
two different behaviors, depending on the relative values of the fluid depth and capillary 
length. At each important step in the calculations, we verify that the shallow water limit 
is recovered. However, the partial differential equations (|2.38|) and (1-14) differ in the order 
of differentiation, with surface tension appearing as a coefficient of the highest derivative 
term in ( |2.38| ); the limit of null surface tension is thus singular. Graphical illustrations of 



the solution are given in Sec. |V], for various values of the dislocation parameter a, and for 
fluid depth larger and smaller than capillary length. An Appendix has some computational 
details. 
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II. WATER WAVES IN INTERACTION WITH A VERTICAL VORTEX 



Equations for an incompressible fluid of equilibrium depth h, free surface h + t](x,y,t) 
with origin of vertical coordinates (z = 0) at the bottom, lying in a (uniform) gravitational 
field g are 

d t V + V-VV = —VP-gz, (2.1) 

P 

V-V = 0, (2.2) 

where V is the fluid velocity, P the pressure and p the (constant) density. 

We neglect viscous dissipation. This is justified if the viscous attenuation time of the 
wave is greater than a typical time for the scattering problem. The attenuation times for 
gravity waves (GW) and capillary waves (CW) of wavelength A are, respectively, f| 



T GW = J»" = (2.3) 



pg 2 ^ T cw = M 2 

2(27r) 4 /ic^' 2(2vr)y 



where \i is the dynamic viscosity of the fluid and the waves phase velocity. In the case 
of water, /i = O.Olg/cms, g = 981cm/s 2 and p = lg/cm 3 . The dispersion relation for 
capillary-gravity waves in a viscous fluid is fairly involved, but the sum of the two times 
gives a good estimate. Numerical estimates for waves of several wavelengths are given in 
Table |I[ The period of the wave is much smaller than the dissipation time in all cases, 
and the travel time on the vortex scale, which is period x (a/A), is also smaller than the 
attenuation time, at least as long as the vortex radius does not become very large. Thus, 
there is a range of values of ka where it is reasonable to neglect dissipation. 

Boundary conditions are that fluid elements at the free surface of the fluid remain there, 
that pressure has a discontinuity that is exactly compensated by surface tension, and that 
there is no vertical velocity at the bottom: 

z = h + rj: V z = d t r] + V ± - V±r) (2.4) 
z = h + T): P = -rV 2 L r/ (2.5) 
z = 0: V z = 0. (2.6) 



with t the surface tension, V± the horizontal velocity and the horizontal gradient. We 
are interested in small perturbations around a steady, axially symmetric, vertical vortex: 

V = U + v v<^U (2.7) 
P = P 0+Pl Pl <P (2.8) 
V = Vo + Vi Vi < Vo- (2.9) 

The vertical vortex is given by the (divergenceless) flow 

U = U (r)6 (2.10) 

in cylindrical coordinates (r,9, z), with (f, z) the unit vectors in the radial, tangential and 
vertical direction respectively. 

We first study the zero order situation, v = 0: 

U ■ VJ7 = — VP — gz. (2.11) 
P 

The 9 component of this equation is an identity. The z component is 

= --d z P o -g (2.12) 
P 



so that 



and the f component is 



Po = -pgz + p {x,y,t), (2.13) 



U 2 1 

= -d rPo . (2.14) 



Given a specific function this is integrated at once. Concerning boundary conditions, the 
third boundary condition ( |2.6|) is satisfied identically. The first boundary condition ( p.4| ) 
says that the surface deformation is independent of polar angle 6, and the second boundary 
condition (|2.5|) gives the free surface t]q in terms of the pressure: 

Po = P9Vo - ^Vi?7o. (2.15) 



Writing v = (u, w) and neglecting terms quadratic in v we have the equations to order 



one: 



(d t + U- V_i_)ti + u ■ V±U = —V±pi 

(d t + U- V ± )w = --d zPl 
V i • u + d?w = 0. 



(2.16) 

(2.17) 
(2.18) 



Similarly, the boundary conditions to order one are || 

z = h + rj: w = (d t + U ■ V_l)^i + u ■ V±rj 
z = h + r): pt = pgr]! - rV^i 
z = 0: w = 0. 



(2.19) 
(2.20) 
(2.21) 



Using (|2.17f) we see that the third boundary condition reads 



z = 0: d z pi = 0. 



(2.22) 



Taking the divergence of Equations (|2.16j ) and ( 2.17 ), and using Q2.18Q gives 



Vipi + d zzPl = -2p(V a U b )(V b u a ). 



(2.23) 



Up to now, these equations are exact for linear surface waves interacting with a static 
vortex. It is the fact that linear waves exist that provides us with another parameter, the 
phase velocity, with which to compare U. We will now simplify the problem by using the 
following two approximations: First, the typical velocity of the vortical flow, Uq, is supposed 
to be much less than the phase velocity of the wave c^, defined in ( |2.4C| ). Second, the 
wavelength A is supposed to be much smaller than a typical length associated with the 
vortex, a. In practice, a will be the core radius of the vortex, and we assume ka ^> 1 where 
k = 27r/A is the wave vector. We will denote formally the small quantities by e. We thus 
assume Uq/c^ = M = 0(e), where M will be called the Mach number in analogy with 
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acoustics, and ka = 0(1/ e), and we search for corrections of order e to the wave equation 
without permanent vortical flow. To get the relative importance of the terms that appear 
in the differential equations, we will use the following estimates: 

V ± / ~-, dth~vfi, V±fi~kf u d,h~kh, (2.24) 
a 

where fo is any scalar quantity refering to the vortical flow, fx any scalar quantity refering 
to the surface waves and v is the wave frequency. We have assumed, as suggested by Eqn. 
(|2.23 ) that length scales for vertical and horizontal variations of surface waves are the same. 



With those estimates, we get from ( |2.17| ) that 

—pi ~ vw =>- p\ ~ pcsw. (2.25) 
P 

Injecting this result in ( |2.23| ), the order of magnitude of the left-hand side is k 2 p\ = k 2 pc^w, 
whereas the right-hand side is of order p(Uo/a)kw = k 2 'pc</,w(Uo/ 'c^)(l/ka); it is thus negli- 
gible, being of order 0(e 2 ), and (|2.23|) is replaced by 

V 2 ±Pl + d zzPl = 0, (2.26) 

which is the same Laplace equation as in the problem of water waves without the vortex; it 
has the big advantage of being autonomous and linear in the pressure so that separation of 
variables can be attempted. 

An estimate of the surface elevation t]q for the vortex flow, may be obtained from ( [2. 15 ) 
and (|2.14|); it reads 



Vo-^tl + ^Y 1 , (2.27) 



9 V a2 



where we have introduced the capillary length l c = yr/pg. For water, r = 74dyn/cm, so 
that l c ~ 0.32cm and the effect of surface tension on the surface deformation of a vortex of 
size a ~ 1 cm, is quite small, of order one per cent. The surface wave elevation, from ( |2.20| ) 



and (|2.25| ), reads 

^-^(l + A^r 1 , (2.2S) 



with k 2 tl of order one. In the following, we take 771 ~ c^w/g, which is numerically inexact 
but adequate for order of magnitude considerations. In ( [2.19| ), the respective orders of 
magnitude of the different terms are 

9 

U-V ±Vl ^^^w = 0(e)(d tVl ), 

9 c<t> 

h. c 2 1 tj2 

u.V ±Vo ^-±±-^ w = 0(e 3 )(d t r ]l ), (2.29) 
and the relevant approximation for fl2.19|) , valid to 0(e), reads 



z = h: w = (dt + U ■ VjjTft. (2.30) 
In the same approximation, we also write 

z = h: px = pgr}! - rV^i (2.31) 
In those equations, we neglect 77 in comparison with h. If we write Pi(h+T]o) = pi(h) + 5pi 



and use fl2.27Q we get 



5pi _ 7]q dpi 
Pi Pi dz 



kr] ~ k— ~ ^f- = 0(e 2 ), (2.32) 



/, 9 cf 



so that Spx is indeed negligible. The same is true for the other boundary equation. 

Let us use now these approximate equations to describe the propagation of surface waves 
in the vortical flow. We will consider almost shallow water waves, that is, the next order in 
the small parameter kh of the calculations of ref. (I). In this limit, the pressure is given as 
a power series in the vertical coordinate z, 

00 

Pl {r,6,z,t) = Y,z n nn{r,0,t). (2.33) 

n=0 

A proper development would consider the dimensionless variable z/L with L a typical length 
scale for horizontal variations. However, ( |2.33| ) is good enough for our purposes. Should 
we wish to have exact results in the deep water limit, kh — > 00, this is the step that would 
break down. Inserting this development in Q2.26D , we obtain the recursion 



S 



which, together with the boundary condition (|2.22| ), gives 



Pi(r,M,t)=£(-l) m ^-^. (2-35) 

m=0 

where we have dropped the index '0' from n for convenience. 

We introduce the notation D t = dt + U ■ Vj_. Taking only the leading order terms in 
the small parameter kh in ( 2.35|) , applying the differential operator D t to (|2.30| ) and taking 



Eqn. ( p.!7|) for z — h, we get 



Dfa = - Vin - ^Vin. (2.36) 
P 6p 



Applying V]_ to ( |2.31| ) and using ( |2.35| ) at the same order, we get 



ghv 2 ±Vl - -vi Vl = -vin - ^vin. (2.37) 

p p 2p 

The surface tension term is considered under the assumption that the capillary length is of 
the same order of magnitude as the depth of the fluid layer. In this case, those two equations 
are valid up to order 0[{kh) 2 ]. It is thus legitimate to replace the pressure by its value at 
order 0(1), II = pgrji, in the term oc Vj_IT, that has the highest derivative. Eliminating the 
pressure in the resulting equations, we get the final result: a dispersive wave equation for 
surface elevation r\\ that is analogous to Eqn. (1-14) of I in the shallow water case. It reads 

ghVlm + (^gh 3 - V A ±m - D 2 tVl = 0. (2.38) 

This equation includes the leading order correction to the shallow water case. It is valid 
under the same assumptions (see ref. [I]) concerning wavelength and fluid velocity. It 
describes the scattering of surface waves over water whose depth is small but not negligible 
with respect to wavelength, when the wavelength is small compared to the vortex size, when 
the velocity of the vortex flow is much less than the phase velocity of the waves, and when 
the waves are of small amplitude. 
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Without the vortex, when U = and d% = D t , plane progressive waves of the form 



77! oc e l(vt - k ^ r ^ 

exist provided frequency v and wave vector k± are related through the dispersion relation 

'rh 1 



V 1 VMS +[i£-lgh'\k\ 



2.39) 



The phase velocity is given by 



c l = ^=9h+ (—- ^gh 3 ] k\ 



and the group velocity c g by 



1 

C4, 



p 



g h+ 2Lj-lgh 3 j k 2 



[2A0) 



(2.41) 



Of course, the dispersion relation ( |2.39| ) is the approximation to order 0[(kh) 3 ] of the well 
known dispersion relation for capillary-gravity waves, 

' rk 3 ^ 



gk + 



tanh kh. 



P 



[2A2) 



The wave dispersion is thus characterized by a dimensionless parameter 8 defined by 

'r h 2 \ 1 



k 



(2.43) 



\pg 3 J 8 

It is positive for h < V3l c , and negative otherwise. We shall consider both cases. In order 
to be consistent with our approximations, namely, that the fourth order term in ( |2.38| ) be 
a small correction to the other two, the absolute value of 8 must be large, and the shallow 
water limit corresponds to \8\ — > 00. For positive 8, the phase and group velocity read 



cl = gh 



1+8 



gh2 + 8 



8 ' s c 8 
whereas for negative values of 8 they read 



(8>0), 



(2.44) 



4 = 9h 



\8\-l 



gh \8\-2 

04, \8\ 



(8<0). 



(2.45) 
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The full dispersion relation ( |2.42 ) for water waves is either convex or concave at small 
depth, depending on the sign of 5. The crossover point h = \f?>l C) derived from the ap- 
proximate relation fl2.39| ), separates two regions of opposite convexity. Both may be exper- 
imentally accessible. The approximation of the hyperbolic tangent is better than 1% for 
kh < 0.5, and better than 5% for kh < 0.8. It is thus easy to stay in the small depth 
limit, tanh(fc/i) m kh — (kh) 3 /3, while keeping the wavelength small in comparison with the 
vortex radius. Using a fluid with high surface tension, like water, leads to a positive 5, that 
is h < \/3l c , whereas the same experiment with a fluid of small surface tension will give a 
negative 6. 



III. SCATTERING OF DISLOCATED WAVES BY A VORTEX 



We will now proceed in close analogy with the calculations of ref. [I]. We consider scat- 
tering of surface waves by a circular uniform vortex with vorticity u and radius a surrounded 
by an irrotational flow. Using polar coordinates (r, 9) , the background flow is given by Q 



U 



\ujrQ if r < a 



if r > a 

iTxr 



(3.1) 



where T = irua 2 is the circulation. Eqn. (|2.38|) will be solved separately for r < a and r > a, 
and the results matched with a continuity condition for 771 and its first three derivatives, 
since ( |2.38| ) is of order four. 

We look for solutions that evolve harmonically (with a single global frequency v) in time, 
and Fourier decompose them in the polar angle 9: 



r/i = Re 



j,(n0—vt) 



(3.2) 



where Re stands for the real part. Inserting this development, and the background flow ( p.l| ) 
for r < a, that is, inside the vortex, in ( ^.38| ), a little calculation shows that the resulting 
equation factorizes exactly as 
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d 2 Id n 2 n , 2 



d 2 Id n 2 „ 



T)ln = 0, 



(3.3) 



with 



{k ± f = L -k 2 8 



-1 ± 



\ 



1 + gh m 1 ' (5>0) ' 



(3.4) 



or 



-fc 2 |<5| 1± 



\ 



4 (i/ - ntu/2) 5 



(5 < 0), 



(3.5) 



The two differential operators in ( |3.3| ) commute, and the four independent solutions of 
this fourth order equation are thus given by the two pairs of solutions of the two correspond- 
ing second order differential equations. 

Taking the shallow water limit 5 — > oo, we get for positive S 



(k, 



gh 



(3.6) 



whereas k_ diverges. Thus k + tends towards the value of k n corresponding to the shallow 
water case (see Eqn. (19) of ref. [I]) as it should, since this case must be recovered as 
a limiting case. The other constant A;_ comes from the fact that ( |3.3| ) is a fourth order 
differential equation, unlike Eqn. (19) of ref. [I]. Its limit for 5 — > oo is singular, reflecting 
the fact that surface tension r multiplies the highest derivative term in the differential 



equation (|2.38| ). The respective role of k + and k_ are exchanged for negative 5. 



From (3.4), we see that when 5 is positive k + is real whereas k_ is imaginary for all 
n, whereas for negative S the two wavevectors k± are real for small n and complex for 
large n. For positive 5, ( |3.3|) has Bessel and Neumann functions as solutions, together 
with hyperbolic Bessel and hyperbolic Neumann functions. The Neumann and hyperbolic 
Neumann functions must be discarded because of regularity at the origin. For negative 5, 
we take Bessel and Neumann functions of a complex argument, and discard the Neumann 
functions to ensure regularity at the origin. Thus 



r/i = Re 



+ b n 



X n (n n r) 
XJK n a) 



(3.7) 
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where the a n and b n are undetermined coefficients in both cases, and where we have intro- 
duced the notation 



h h — iK X — T 

^ny rv — — m nj yi -n — ± ni 



(6 > 0), 



(3.8) 



k n , k+ — K n , X n — J n , (5 < 0). 



(3.9) 



Outside the vortex, for r > a, using U = T/(2irr)9 and the decomposition ( |3.2| ), and 
dropping terms of order M 2 , we get that ( |2.38| ) reads 



Vln = 0, 



where, using the dispersion relation (|2.39| ), we define the two constants 



A=(5 + l)k\ B = — 



Sk 2 2Tvn 



gh 2n 



and the linear differential operator 



r = — - — 

dr 2 r dr 



(3.10) 



(3.11) 



(3.12) 



This equation may be written in the factorized form 



o + o^ ln = 0, 



(3.13) 



with 



mi n 
± = C-^ + q 2 ± 



(3.14) 



provided the unknown coefficients m+, m_, q + and g_ satisfy the following relations: 



r 
r 2 ) 



r 3 dr 
1 



q 2 + + q 2 _ = -5k 2 
qW- = -A, 



m + q_ + m_q + = —ok n — n, 
m 2 , + m 2 _ = 2n 2 , 



m = n , 



OO 0/1 

m,m_ — 4m_ = n —An 



(3.15) 
(3.16) 
(3.17) 
(3.18) 
(3.19) 
(3.20) 
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Here we have indicated on the left the portion of the linear differential operator that leads to 
each condition. There are six equations for only four unknowns and obviously they cannot 
be simultaneously satisfied in general. The last two equations, ( |3.19| ) and ( |3.20| ), correspond 
to terms that are negligible at large distance from the vortex. If we compare them to C 2 , 
they are smaller than 1//3 2 because r > a. Accordingly, we are justified in ignoring these 
two equations in our approximation (5 ^> 1 and we solve (|3.15| - ^T8|) , which gives 



q ^ = k 2 >0, q 2 _ = (iqf = -k 2 {\ +5) < 0, (m ± ) 2 = n 2 ± 2na, (5 > 0), (3.21) 

q\ = q 2 = (|,J| -l)k 2 > 0, q 2 _ = k 2 >0, (m±) 2 = n 2 T 2na, (5<0), (3.22) 

_Yu 1 _ IV 1 

a= 2^gh + 2(r/p - gh 2 /3)hk 2 ~ 2^ 7^ (3 ' 23) 

where we used Q2.41p to write the last equality. It is important to note that the index 
n 2 + 2na is always associated to the incident wavevector k. We will comment further on 
this result after Eqn. ( p.27|) . From now on, we define 



Vn 2 + 2na, m_ = Vn 2 — 2na, (3.24) 



so that in the negative S case m_ (resp. m + ) is associated with q + (resp. g_), as shown by 
O). 



The dimensionless parameter a is defined in close analogy with I. We can write a = 
M/3(c^/cg), which may be of order 1, whith M < 1 and j9 > 1. This parameter has the 
same physical interpretation as in the shallow water case (see below) : it gives the amount 
of dislocation for the wavefronts far from the vortex. This calculation provides an explicit 
confirmation, in a perturbation expansion near the shallow water case, of the intuitive result 
of Berry et al. @ . 

The two differential operators 0± in ( |3.13| ) do not commute. Using the usual notation 



for the commutator of two operators, we get 



[o + , oa = k - ™ 2 -) k M = 4 K - --) < 325 ' 
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which is small, of the same order as the neglected terms fl3.19| , |3.20|) , and will also be 



neglected. Thus, in the same approximation, for positive 5 the solution of ( |3.13| ) is a linear 



combination of Bessel, Neumann, hyperbolic Bessel and hyperbolic Neumann functions, 
because q + is real and g_ is imaginary. Since the hyperbolic Bessel function diverges at 
infinity, it must be discarded. For negative 5, the solution is a linear combination of Bessel 
and Neumann functions of argument kr and qr. Since the wave number q is that of a 
scattered wave, we discard the Bessel function of qr, keeping only the outgoing wave from 
the vortex. 

Following Berry et al. || as in [I], we write the surface elevation outside the vortex in 
the form 

i]i = He(i]AB + t]r), (3.26) 

where 

„ _ „ Jm+{kr) u n8 -ut) ( o 97 n 

n J m + \P) 

with (3 = ka. It does not depend on the sign of 8, which is physically obvious because 
the amount of dislocated wavefront is linked to the circulation of the vortex, not to the 



curvature of the dispersion relation. Thus m + = + 2na is always the index of the 
functions involving the wavevector k. The other term of (|3.26| ) depends on the sign of S, 



which is also physically clear since they represent the wave scattered by the vortex. They 
read 

Depending on the sign of S, we have the following definitions : 

Y m _=K m _, q = kVTT5, (6>0) (3.29) 



Y m _=Hi_, q = kJ\5\-l, (6<0) (3.30) 
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The coefficients a n , b n , c n , d n and e n are defined so that they denote the amplitude of the 
wave components at the vortex boundary r = a. In order to obtain these coefficients, and 
since the equation ( |2.38| ) is of order four, the continuity of fji and its first three derivatives 



at r = a is required. This gives four relations: 

Q"n 4" b n d n e n = c n , (3.31) 

Ji<{k n a) 2 KM _ ^'{(3) _ Y»_(ga) _ J^(fl 

ankn J n {k n a) +bnKn X n M dnk m m+ (f3) Y m _(qa) Jm+ifiY [ ' 



where X (resp. K) is defined either by Q or Q [resp. (ggg) or (|3~30|) 1. 



The fifth, and last, boundary condition comes from the asymptotic behaviour of rj at 
infinity. We require that the asymptotics of tjab coincides with the dislocated wave incident 
from the right plus outgoing waves only. Exactly in the same way as in I, this leads to 

-t) m+ , (3.35) 



It is important, in order to use this result, that either the coefficient q+, for positive 5, or 
gi, for negative 5, in ( |3.13| ) be equal to k 2 , and that they be associated in each case to m+. 



Otherwise it would have been impossible to recover the dislocated wave, which is a crucial 
physical requirement for the solution because we need q = k to be a possible result of the 
factorization. This fact fully justifies the factorization in (|3.13|) . 

The solution of system (|3.31| - j3l34D is thus known in principle, but it is not too illumi- 
nating, and it will not be displayed here. Details of the calculations may be found in the 
Appendix. To illustrate the solution, we use Mathematica |9j in order to do the calculations 



that will be indicated in Sec. |V|. Thus it is sufficient to have the coefficients expressed 
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as ratio of 4 x 4 determinants, as in ( |A9|) . We insist on the fact that the solution may be 
innacurate at a few wavelengths away from the vortex, because of the approximate character 
of the factorization ( |3.13| ) . 

Let us discuss the asymptotic behavior of the solution ( |3.27| , |3.28| ) for r — > oo. The case of 
t]ab is completely similar to the shallow water case. Indeed, the index of the Bessel function, 
m + {n) = \Jn 2 + (Const.) x n, has exactly the same structure as m{n) in [I]. An important 
consequence is that the parameter a = /3M(c^/c g ) has the same physical significance as in 
the shallow water (or acoustics) case : it quantifies the dislocation of the wavefronts in the 
forward direction, at large distances from the vortex. Other results may also be transposed 
in a straightforward fashion, and the asymptotics of 7]ab i r ,@) for large r is still given by 
Eqn. (38) of [I], with the proviso that the function G{6, — ir/2) be replaced by 



Gdw 

where m + is given by (|3.24j) and 



,t)= E ^ 

\n\<N 



7r/2) _ e im i d (i-7r/2) 



(3.36) 



Tu 1 
n H 



(3.37) 



27T C<pC g 

The asymptotics of t]r depends on the sign of 5. If 5 is positive, the hyperbolic Bessel 
function does not contribute to the scattered far field, because |I| K p (z) ~ e~ z / ' y/z for 
large z. Thus 



1/2 



t i(kl — vi) 



d 



( 'Em " 



n i{n9— 7rm + /2) 



(3.38) 



Mkr) n H L+(^ 

In the next section, we will compare the correction to the Aharonov-Bohm scattering am- 
plitude in the case of shallow water waves, given by Eqn. (43) of [I], and the correction for 
dispersive water waves, which reads 



dr, 



fuwiO) = G DW (6, -vr/2) +2j2--^-e M (-i) 

n h mAP) 



m_|_ 



(3.39) 



If S is negative, we must take into account the two outgoing Hankel functions. We get 



VR 



nikr 



1/2 



E 



d -inm+/2 



+ 



<3 r , C 



-i7rm_/2 



And 



(3.40) 
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and the correction for dispersive water waves now reads 



/w(0)=GW(0,-7r/2)+2X: 



d n (-i) m + e n (-i) 




(3.41) 



n 



IV. NUMERICAL EXAMPLES 



The solution to the scattering problem of surface waves by a uniform vertical vortex 
depends on four dimensionless parameters. A first set includes the dimensionless vortex 
radius (3 3> 1 and the dislocation parameter a = /3M(c^/c g ) which quantifies the wavefront 
dislocation. They already appeared in [I], with the same definitions and physical interpre- 
tations. A third parameter is the dimensionless capillary length i = kl c , and the last one is 
the dimensionless depth kh. In order to simplify somewhat the discussion, we use the single 
dimensionless parameter 8, defined in ( |2.43| ), in place of the depth and capillary length. As 
an example, we take 6 = 8, which may correspond, for example, to h = l c , kh = a/3/4, and 
5 = —8, which may correspond to h = 3l c , kh = 3/4. In both cases, the hyperbolic tangent 
in (|2.42|) is approximated to better than five percent by the two leading terms, the ones we 
are keeping, in its series expansion. 

Scaling radial distance with the vortex radius, r' = r/a, the analytical expression of the 
surface displacement is summarized as follows, depending on the sign of 5. Inside the vortex 
(0 < r' < 1) we have 



771 = Re r] c . 



For positive values of 5, 




(4.1) 



where we have defined the following dimensionless wave numbers: 




2 



1/2 



1 — n 



(5 2 l + 5 



a 2 + 5 



) 



(4.2) 
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n n a = (p n = (3 d- 



1 + 



, 1 + 5/ a 2 + 5' 



1/2 



For negative values of 5, 



whith new dimensionless wave numbers: 



/c n a = (f) n = \ —f3 



( 



V 



1/2 



<5 2 



K n a — (p n — 



1 + 



V 



1/2 



Outside the vortex (r' > 1) 



where, whatever the sign of 5, 



?7i = Re(?7AB + 77^), 



^B = E(-0 m+ ^ + (/3r')e 



l\ i(n6—vt) 



For positive values of 5, 



1 < ' K m _ (j3VT+5ry 



where we have used the fact that 



qa = PV1 + 5. 



For negative values of 5, 



HUPr') , H^lpJ\6\-lr> 



VR = T, K m (os + e 



j{n6—vt) 



where we have used the fact that 



qa = 



I. 
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Some details of the calculations of the coefficients a n ,b n , d n and c n are given in the Appendix. 
As an illustration, absolute values of the coefficients a n and b n (resp. d n and e n ) are plotted 
in Fig. 1 (resp. in Fig. 2), for positive 5, and in Figs. 3 and 4 for negative 5. 

Since convergence of the series expansions for tjab and tjr is not uniform, the number of 
terms to keep in the numerical evaluation of the infinite series depends on the value of r' . 
In practice, the convergence of the series is comparable to the case of ref. [I], and we use 
roughly the same number of terms. We compute the patterns of the surface displacement 
in the region \y'\ < 5[(x',y') = (r' cos 8 , r' sin 8)} by the finite sum of (|4.1| , |4.4j ) and 



1.10) with \n\ < 50 for (3 = 10 and |n| < 30 for (3 = 5, but we keep more terms, \n\ < 90 in 



Let us first consider the case of positive 5. Fig. 5 shows the resulting displacements 
for 5 = +8, (3 = 5 and a = 0.5, 1, 1.5, 2, and Fig. 6 for the same values of a and 5, but 
(3 = 10. The dislocation of the incident wavefronts by an amount equal to a is clearly 
visible. The outward travelling scattered wave is also visible. The interference patterns 
between scattered and incident wave is very similar to the corresponding pictures of [I], for 
the same values of (3 and a. This is confirmed by the comparison of the scattering cross 
section displayed in Fig. 9, as discussed below. Taking into account the dispersion greatly 
modifies the numerical value of a, but other corrections are small in the case of positive 5. 

In the case of negative 5, the interference pattern is strongly modified. This is illustrated 
in Fig. 7, where we plot the surface displacement for 5 = —8, (3 = 5 and a = 0.5, 1, 1.5, 2. 
The spiral wave which is clearly seen for a = 2 (Fig. 7-d) is very different from the corre- 
sponding figure of [I] (see Fig. 2-d of [0). For larger values of (3, shown in Fig. 8 for which 
(3 = 10, whith same values for S and a as before, the pictures are much more similar to the 
shallow water case. 

This is confirmed by the plot of the absolute value of the correction to the Aharonov- 
Bohm scattering amplitude in both cases. The dashed line in Fig. 9 shows this correction in 
the shallow water (non dispersive) case, and the solid line shows the same correction in the 
dispersive case, for a positive 5 = +8. Both corrections are almost the same, in agreement 
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with the results of Fig. 5 & 6. The same is done for a negative 5 = —8 in Fig. 10. The graphs, 
this time, differ enormously. In the dispersive case, the scattering is much more isotropic, and 
very different in amplitude. An obvious, but somewhat formal, explanation of this difference 



is the supplementary function H] n _ (PJ\5\ — 1) in ( |3.41| ) compared to ( |3.39| ). Also, the index 



of this function is m-(n) which takes imaginary values for small positive n, rather than for 
small negative n as m + (n). This implies that the partial amplitudes for exp(— in8) and 
exp(m#) are much more similar to each other than in the shallow water case, for which 
m_ is absent, and also more similar than in the case of positive 5, where the decrease of 
the corresponding function is exponential. The appearance of an algebraically (oc 1/y/r) 
decreasing amplitude associated to m_ in the negative 5 case restores the symmetry of the 
scattered wave. 

A more physical explanation is as follows: Consider a plane wave incident from the right 
on a counterclockwise vortex, as in figures 5-8. Above the vortex, the wavefront velocity 
is increased by advection, whereas it is decreased below the vortex. Consequently the wave 
fronts should bend towards the bottom of the picture, as they do. The other effect of the 
vortical flow is to add a wavelength below the vortex, which means that the wavenumber 
k decreases below the vortex. Since kh is assumed to be small, for positive 5 the phase 
velocity increases with k. The phase velocity is thus smaller below the vortex than above, 
which enhances the bending of the wavefronts, and reinforces the strong asymmetry in the 
interference pattern of Fig. 5 & 6, and in the scattering amplitude of Fig. 9. For negative 
S, the phase velocity decreases with k, and the effect of the dislocation is to make the phase 
velocity smaller for the part of the wavefront above the vortex. This effect balances the 
effect of advection, and we understand why the interference pattern (see Fig. 7 & 8) and 
the scattering amplitude (see Fig. 10) are much more symmetric than in the positive 5 
case, or than in the nondispersive (5 = 0) case. It is also reasonable that this effect should 
be more important for (3 = 5 than for (3 = 10, because the relative variation in k due to 
the dislocation is greater in the former case. The spiral waves are observed for negative 
S because the interference pattern almost keeps rotational symmetry while smoothing the 
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wavefront dislocation in the forward direction. 

As a last illustration, we compare the wave patterns predicted by the shallow water 
approximation, and by its first correction in powers of the fluid depth, in an experimentally 
accessible situation. We suppose the fluid to be pure water, of depth 1 mm.; the vortex 
radius is 1 cm., and the wavelength is 2 cm. Thus kh = vr/10, and the approximation of 
the dispersion relation is excellent. The price to pay is the rather small value (3 — 7r. We 
take the vortex circulation such that a = 1 in the shallow water approximation, for which 
c = \fgh = 9.9 cm/s. Taking the dispersion into account, we get 5 = 1.4, = 13.0 cm/s, 
c g = 18.4 cm/s and a = 0.41. The difference in the respective numerical values of a is the 
leading effect. The result is shown in Fig. 11. To obtain quantitative agreement with an 
experimental situation, it may be sufficient to keep the shallow water approximation, but 
with the actual value of a obtained in the dispersive case ( p.23|) . It should be interesting to 
use a fluid with a very small surface tension, in order to obtain a negative value of 5 while 
keeping a small value of kh, for which the wave pattern should be extremely different from 
the shallow fluid layer approximation. All calculations were performed using Mathematica 



V. CONCLUDING REMARKS 

We have computed the surface displacement due to a dispersive surface wave interacting 
with a vertical vortex when the vortex core performs solid body rotation, the wavelength 
is small compared to the vortex core radius and the particle velocities associated with the 
wave are small compared to the particle velocities associated with the vortex. When the 
parameter a = uT /27rc^Cg is of order one or bigger, the wavefronts become dislocated. This 
parameter depends, in the dispersive case, both on the phase and group velocity of the wave, 
and tends smoothly toward the result of [I] in the nondispersive limit. We thus give a proof, 
in a perturbative fashion, of the heuristic derivation of Berry et al. |2j . 

Formally, we proceed perturbatively around the shallow water limit, to obtain a fourth 
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order partial differential equation for the surface elevation associated with the surface wave. 
However, apart from some technical details, the solution is very similar to the nondispersive 
case. The scattered waves interact strongly with the dislocated wavefronts and produce 
interference patterns. A dimensionless parameter 5 quantifies the relation between fluid 
layer depth h and capillary length l c . When it is positive (h < y/3l c ) the wave pattern is 
similar to the shallow water case. When it is negative, for large values of the circulation, 
the wave pattern is very different. 

We hope that the calculations in the dispersive case will help the comparison with exper- 
iments. Our calculations are valid when the approximation tanhkh ~ kh — (kh) 3 /3 holds. 
This is a restrictive condition, but we believe that, once dispersion is correctly taken into 
account in the definition of a, the wave pattern given by the nondispersive case should be 
roughly similar to the observations. Some discrepancies are expected when the fluid depth 
is greater than the capillary length (h > v^3Z c ), but in this case it is necessary to correctly 
approximate the hyperbolic tangent by the first two terms of the series, and in practice the 
effect should be observable for a fluid of small surface tension (thus small capillary length) 
only. 
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APPENDIX A: COMPUTATIONAL DETAILS 



In this Appendix, we explain how to calculate the coefficients in the series representations 
of Sec. [TTT]. We also discuss briefly the convergence of the series involved in our solution 
( j3.7| ) and ( |3.28|) , since the sum in Q3.27| ) obviously converges as in the shallow water case [I] . 

To get the solution of eqns. (|3.31| - |3~i?4T) , we introduce the column vectors V x of the 
coefficients of the unknown x n . They involve multiple derivatives of the Bessel functions, 
and may be simplified with the help of the following formulae \TU\: 

i , „ . r _ _ 2 P 



j p-i 



j p+i 



+ 



< I 



K' 



■(K p . 1 + K p+1 ), 



p-i 



I + 2P I 
1 p+i "I 1 pj 

z 



(Al) 



= K, 



2p 



p+i 



K„ 



V 



where z is the argument of the function, p its index, and Z represents either J or H l . 
The column vectors thus read 



/ 



V a = 



n - (j) r , 



1 

J, 



n+lv 



J 1 1 



n(n -1) - 4>l + (fi n 



J, 



J n ((fi n ) 



V („ _ i)[ n (n - 2) - eg] + <f> n [<f>l ~ + 2)] 



J, 



n+lWn) 
Jr,.(d>r,) 



(A2) 



in the negative S case, and the same expression with (fi n in the positive 5 case, 



V h = 



1 



\ 



n(n -l) + ip 2 n -ip n 
V (n - 1) [n(n - 2) + Vn ] + ^ + n 2 + 2) + 1 



(A3) 



in the positive 5 case, or 
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Vh 



1 

~ Jn+l{<fn) 

n-ip n ——— 



n(n -l)-(p n + <f> 



~ Jn+l(fn) 



Jn(<Pr. 



V („ _ i)[ n ( n _ 2 ) - + £ n [$ _ ( n 2 + 2)] 



Jn{<fn) 



in the negative 5 case, 



V A 



— m+ + /3 



#m + +lC#) 



. m+(m+ _ 1)+/32 _,M 

V (m+ - 1) [/3 2 - m + (m+ - 2)] + /?(w 2 + 2 - /3 2 ) 



-771. 



+ PVTT5 



-l 
a: 



m_ +1 



K„ 



-m_(m_ - 1) - /3 2 (1 + 5) - /3^m 



1 + 5 



1 + 5 



1 + 5 



(1 - m_) /3 2 (1 + 5) + m_(m_ - 2) 
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V 



in the positive 5 case, or 
/ 

- m _+(3^\5\ 
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1 + 5 



V e 
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rrl 
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(A4) 



(A5) 



(A6) 



(A7) 



in the negative 5 case. The column vector V c for the coefficients of c n in the right hand 
member of ([OT]-[03]) reads 
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-i) 



m_|_ 



/ Jm+iP) \ 

m+Jm+{(3) - /3J m++ i((3) 
[m + (m + - l)-P 2 ]J m+ (P)+pj m++ i(P) 
(m+ - l)[m+(m+ - 2) - (3 2 }J m+ {(3) + 

V +P[p 2 -(m 2 + + 2)]J m++1 (P) j 

Now, formally, the coefficient a n (say) is expressed as the ratio of two determinants, 

\v c v b v d v e \ 



(A8) 



(A9) 



\v a v b v d v e \ 

and is a function of n, a, f3 and 5. 

Let us discuss briefly the asymptotic behavior of those coefficients. For large n, m + ~ n, 
and we can safely assume that Z n+1 /Z n = 0(1), in the case of a constant argument, /3, 
(3\/l + 5, (3\J\5\ — 1, or in the case of an argument of order n, such as n or ip n . This is 
demonstrated for ordinary Bessel functions in [I], and may be deduced in the same fashion 
for hyperbolic Bessel functions, which have similar behaviors for large values of the index 
| TTJ . From the calculations of ref . [I] , we easily get the dominent behavior for large n of the 
vectors V x , and thus of the coefficients x n . If x n is not equal to c n , 

( 1 \ 



V, 



n 



\n 3 ) 



(A10) 



so that the determinant in the denominator of any expression such as (|A9|) behaves like 



\v a v b v d v e 



~n 6 . 



(All) 



We also have that 



( 1 ^ 




( 1 ^ 


n 
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UP) ~ 




n 2 




n 2 


\n 3 J 




\n 3 J 



n \n 



(A12) 
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so that the determinant in the numerator of (|A9f) behaves like 



\V c V b V d V e \ ~ n 5 / 2 f£) , (A13) 

and the coefficients x n decrease sufficiently quickly to ensure the convergence of the series in 
( P-7| ) and ( |3.28| ). The convergence is uniform for r < a, because the support of the functions 



is compact, but not for r > a. 
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FIGURES 

FIG. 1. Absolute magnitude of coefficients a n (a) and c n (b) versus n, in a log-linear plot, 
for a positive 5 = +8, and for values of dislocation parameter a and dimensionless radius (3 
(a, (3) = (0.5, 10), denoted by dots, (a, (3) = (1.5, 10), denoted by empty circles o, (a, (3) = (0.5,5), 
denoted by filled circles • and (a, [3) = (1.5, 5), denoted by empty squares □. Note the asymmetry 
with respect to n — > — n 

FIG. 2. Same as figure 1, for the coefficients d n (a) and e n (b). 

FIG. 3. Same as figure 1, for a negative 5 = —8. 

FIG. 4. Same as figure 2, for a negative 5 = —8. 

FIG. 5. Density plot of the surface elevation for the total wave patterns for 5 = +8 and (3 = 5, 
for several values of a. Respectively a = 0.5 : (a), a = 1 : (b), a = 1.5 : (c), a = 2 : (d). The 
greyscale is linear with surface amplitude (arbitrary units). The dark ring indicates the vortex 
location, and the vortex rotates counterclockwise. The wave is incident from the right. The box 
size is 10 x 10 in units of a, the vortex radius. 

FIG. 6. Same as figure 5, for 5 = +8, (3 = 10, and several values of a. Respectively a = 0.5 : 
(a), a = 1 : (b), a = 1.5 : (c), a = 2 : (d). 

FIG. 7. Same as figure 5, for S = —8, (3 = 5, and several values of a. Respectively a = 0.5 : 
(a), a = 1 : (b), a = 1.5 : (c), a = 2 : (d). 

FIG. 8. Same as figure 5, for 5 = —8, (3 = 10, and several values of a. Respectively a = 0.5 : 
(a), a = 1 : (b), a = 1.5 : (c), a = 2 : (d). 
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FIG. 9. Polar plot of the absolute value of the correction to the Aharonov-Bohm (i.e. point) 
scattering amplitude, in the case of non dispersive waves (dashed line) and in the case of dispersive 
waves (solid line), in the case of a positive 5 = +8, for (5 = 5 and a = 0.25 : (a), a = 0.5 : (b), 
a = 1.0 : (c), a = 1.25 : (d). The vortex location is marked by the large dot; the vortex rotates 
count er clo ckwise . 

FIG. 10. Same as figure 9, but for a negative 5 = —8, for j3 = 5 and a = 0.25 : (a), a = 0.5 : 
(b), a = 1.0 : (c), a = 1.25 : (d). The vortex location is marked by the large dot; the vortex rotates 
counterclockwise. In this case the dispersive wave scatters very differently from the nondispersive 
wave. 

FIG. 11. Density plot of the surface elevation for the total wave pattern calculated in the 
shallow water approximation, (a), and to first order in fluids depth (b). The greyscale is linear 
with surface amplitude (arbitrary units). The dark ring indicates the vortex location, and the 
vortex rotates counterclockwise. The incident wave comes from the right. The box size is 20 x 20 
in units of a, the vortex radius. j3 = -it in both cases, but a = 1 in the shallow water case (a), and 
a = 0.41, 5 = 1.4 in the dispersive case (b). 
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TABLES 





C-G waves 


S-W waves 


A (cm) 


0.1 


0.5 


1 


2 


c<£ (cm/s) 


68 


32 


25 


30 




10- 7 


0.0018 


0.08 


0.6 


T cw (s) 


0.013 


0.3 


1.26 




Wave period (s) 


0.0015 


0.0156 


0.04 


0.07 



TABLE I. Attenuation times for capillary-gravity (C-G) waves and shallow water (S-W) waves, 



compared to the period of the wave. The attenuation times are defined by (2.3) 
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Coste et al, Figure 
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C. Coste et ai, Figure 
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C. Coste et ai, Figure || 
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C. Coste et ai, Figure [8] 
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